The general Two-Higgs doublet extensions of the SM: a saucerful of secrets 
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I. INTRODUCTION 
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We discuss the most general formulation of the Two-Higgs doublet model, which incorporates 
flavor changing neutral scalar interactions (FCNSI) and CP violation (CPV) from several sources. 
CP violation can arise either from Yukawa terms or from the Higgs potential, be it explicit or 
spontaneous. We show how the model, which is denoted as 2HDM-X, reduces to some versions 
known in the literature (2HDM-I,II,III), as well as some of their variants (top, lepton, dark) denoted 
here as 2HDM-IV. We also discuss another limit that includes CPV and Yukawa four textures to 
control FCNSI, which we denote as 2HDM-V. We evaluate the CPV asymmetry for the decay 
h — > bcW , which may allow to test the patterns of FCNSI and CPV, that arises in these models. 

PACS numbers: 12.60.Fr, 12. 15. Mm, 14.80.Cp 
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Despite the success of the Standard Model (SM) in the gauge and fermion sectors, the Higgs sector remains the 
least tested aspect of the model [l| , which leaves the puzzles associated with the mechanism of electroweak symmetry 
breaking (EWSB) still unsolved. On one hand, the analysis of raditive corrections within the SM points towards 

II ' the existence of a Higgs boson, with a mass of the order of the EW scale, which in turn could be detected at the 
Qh LHC [6, 7]. On the other hand, the SM is often considered as an effective theory, valid up to an energy scale of 
0(TeV), that eventually will be replaced by a more fundamental theory Q, which will explain, among other things, 
the physics behind EWSB and perhaps even the origin of flavor. Many examples of candidate theories, which range 
from supersymmetry [9Hll| to strongly interacting models [l2l . [l3j as well as some extra dimensional scenarios [T3 - 
llql . include a multi-scalar Higgs sector. In particular, models with two scalar doublets have been studied extensively 
[17h19| . as they include a rich structure with interesting phenomenology [20l - |22j . 

Several versions of the 2HDM have been studied in the literature [23| . Some models (known as 2HDM-I and 2HDM- 
II) involve natural flavor conservation (24|, while other models (known as 2HDM-III) (23|, allow for the presence of 
flavor changing scalar interactions (FCNSI) at a level consistent with low-energy constraints [25[ . There are also some 
variants (known as top, lepton, neutrino), where one Higgs doublet couples predominantly to one type of fermion 
[27| . while in other models it is even possible to identify a candidate for dark matter (26|. The definition of all these 
models, depends on the Yukawa structure and symmetries of the Higgs sector [281 - 433 ] , whose origin is still not known. 
The possible appearance of new sources of CP violation is another characteristic of these models [33[ . 

In this paper we aim to discuss the most general version of the Two-Higgs doublet model (2HDM), which incorpo- 
rates flavor or CP violation from all possible sources [3J-|36(. We also discuss how the general model, denoted here 
as 2HDM-X, reduces in certain limits to the versions known as 2HDM-I,II,III, as well as some of their variants which 
we shall name as 2HDM-IV and 2HDM-V. The logic of the naming scheme that we adopt here, consists in identifying 
distinctive physical characteristics that can be associated with the models and have sufficient merit to single out them. 

Within model I (2HDM-I) where only one Higgs doublet generates all gauge and fermion masses [6|, while the 
second doublet only knows about this through mixing, and thus the Higgs phenomenology will share some similarities 
with the SM, although the SM Higgs couplings will now be shared among the neutral scalar spectrum. The presence of 
a charged Higgs boson is clearly the signal beyond the SM. Within 2DHM-II one also has natural flavor conservation 
[2~4j . and its phenomenology will be similar to the 2HDM-I, although in this case the SM couplings are shared not 
only because of mixing, but also because of the Yukawa structure. On the other hand, the distinctive characteristic 
of 2HDM-III is the presence of FCNSI, which require a certain mechanism in order to suppress them, for instance 



one can imposes a certain texture for the Yukawa couplings 
couplings [3|. Within all those models (2HDM I,II,III) [3! 



371. which will then predict a pattern of FCNSI Higgs 
4l| . the Higgs doublets couple, in principle, with all 



fermion families, with a strength proportional to the fermion masses, modulo other parameters. 
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There are also other models where the Higgs doublets couple non-universally to the fermion families, which have 
also been discussed in the literature [13, IH, HJ we shall denote this class of family non-universal models as 2HDM-IV. 
In principle, the general model includes CPV, which could arise from the same CPV phase that appears in the CKM 
matrix, as in the SM, from some other extra phase coming from the Yukawa sector or from the Higgs potential [44 |. 
However, in order to discuss which type of CP violation can appear in each case, we shall use the label 2HDM-V to 
denote the class of models which, besides containing a generic pattern of FCNSI, modulcd by certain texture, will 
include new sources of CPV as well. 

Our formulation of the 2HDM-X, which is discussed in detail in section 2, relies on the Higgs mass-eigenstates basis 



45j. It seems to us that this is more appropriate in order to relate the low-energy constraints on the parameters of 



the models, with the predicted high energy signatures to be searched at future colliders. The different models can be 
characterized by a set of invariants, signaling the possible appearence of CP violation [33], either from the bosonic or 
fermionic sectors. Section 3 contains the discussion of general Yukawa couplings, with spontaneous or explicit CPV 
parameters for the Higgs bosons. We then show in section 4, our evaluation of the CPV asymmetry for the decay 
h — » bcW, which may allow to test these patterns of FCNSI and CPV at a future colliders. Finally our conclusions 
are included in section 5, while some technical details are included in the appendices. 



II. A GENERAL FORMULATION OF THE 2HDM AND ITS LIMITING CASES. 



The Two-Higgs doublet extension of the SM includes two scalar doublets of equal hypercharge, denoted by: $1.2 = 
(4>i2i ( l ) i 2) 1 Depending on the Yukawa matrices Y^ 2 (<7 = u : d) that are allowed, one defines the particular versions 
of the 2HDM. FCNSI appear at tree level when more than one Higgs doublet couples to both types of quarks (it and 
d) 0, and a certain mechanism should be invoked in order to bring them under control. The CP properties of the 
Higgs boson depend on the symmetries of the potential [44[ • 

In order to clarify the discussion of the many models that have been presented in the literature, we shall present 
in the next subsections, a classification scheme for these models, which have different patterns of FCNSI and CP 
properties. We shall discuss in this paper first the most general formulation of the 2HDM-X, and will consider 
the 2HDM versions usually discussed in the literature, which are known as 2HDM-I,II,III, as well as some variants 
(2HDM-IV). We shall discuss then in detail, some cases that have not been considered before, which we denote as 
2HDM-V. Although the 2HDM-X suffers from the FCNSI problem, we shall discuss it first in general terms, without 
referring the specific mechanism that is used to address the problem, which will be done later in this section. 



A. A classification of models 



We shall define here the different types of models according to their Yukawa structure, the Hermiticity of the 
Yukawa matrices and the CP properties of the bosonic Higgs sector. Thus, the most general version of the 2HDM is 
defined through the following assumptions: 

i) In principle we allow each Higgs doublet to couple to both type of fermions, expecting that some particular 
structure of the Yukawa matrices is responsible for the suppression of flavor changing neutral scalar interactions 
(FCNSI). 

ii) The Yukawa matrices are allowed in general to be non-Hermitian, i. e., Yfi 7^ Y^ (/ =u, d, 1, i = 1, 2). The limit 
when the Yukawa matrices are hermitic defines a particular version of the models. 

iii) The Higgs potential admits in principle both spontaneous or explicit CPV. 

Then the known limiting models (2HDM I, II, III), arc obtained by relaxing some of those assumptions, namely: 

1 . The 2HDM-I [l?], HE S3 is defined by considering that only one Higgs doublet generates the masses of all types 
of fermions, as it happens in the SM. This type of model can be obtained by assuming an additional Z2 discrete 
symmetry. Under a variant of this model, where the second doublet does not mix with the first doublet, it is 
possible to identify a neutral scalar as a dark matter candidate (23|, which makes it very attractive. 

2. For the so called 2HDM-II [13, E3, IHl , each Higgs doublet couples only to one type of quark, and then FCNSI 
do not appear at tree level. A variant of the Z2 discrete symmetry is considered here, similarly to the case of 
the 2HDM-I. Two limiting cases can be considered, namely: the 2HDM-IIa, with CP conserving Higgs sector, 
and THDM-IIb where the Higgs sector ir CP violating [49] . This model is also attractive because it corresponds 
to the Higgs sector of the Minimal supersymmetric standard model (MSSM), at tree- level [g. 
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Model type 


Up quarks 


Down quarks 


Charged leptons 


Neutral leptons 


2HDM-I 


#1 


Hi 


#i 


#i 


2HDM-II 


#2 


m 


m 


#2 


2HDM-III 


#1,2 


#1,2 


#1,2 


#1,2 


2HDM-IV 


#1 


#1 


#1 


#2 



TABLE I: Higgs interaction with fermions for 2HDM types. 



Model type 


FCNC 


Hermiticity 


Higgs sector CP 


I 


X 






II 


X 






III 


/ 


/ 


/ (CKM) 


Va 


/ 


/ 


X 


Vb 


/ 


X 


/ 



TABLE II: Symmetries under the different types of 2HDM's 



3. Within the 2HDM-III [38[, one considers all possible couplings among the Higgs doublets and fermions in the 
Yukawa sector; thus, it is possible to have FCNSI in this case. According to the extended classification that we 
try to motivate here, we shall also assume that within 2HDM-III the Yukawa matrices are Hermitic, whereas the 
Higgs potential is CP conserving. Thus, CP violation only arise from the CKM phase. A particular version of 
this model, widely studied in the literature, assumes Hermitic Yukawa matrices with 4-textures [371 ] . which has 
under control the FCNSI problem (38|. It also happens that when one considers loop effects within the MSSM, 
its Higgs sector also becomes of type III which again makes attractive this version of the 2HDM. Although 
it is not often explicitly stated, we shall consider that within 2HDM-III the Higgs doublets couple in principle 
with all three families of quarks and leptons. 

4. Family non-universal assignments are also possible, for instance we can have models where one doublet couples 
to all types of quarks, but a second doublet only couples to the 3rd family (27[. Several possibilities have been 
considered in the literature, which we denote as 2HDM-IV, de pen ding on whether the second doublet couples 
to the whole third family or only to the top (2HDM-IV-t) [42|, |43| . We shall also include within this category, 
those models where one doublet couples only to charged leptons or to neutrinos 0, |43| . 

The properties of these models are summarized in tables HI [II] Table [TT] shows the assumptions for the different 
types of the 2HDM which are considered in this work. Besides the cases, that have been discussed in the literature, 
we can define still another class of models, which we denote as 2HDM-V, where besides having FCNSI at tree level, 
also include extra sources of CP violation, either from the Yukawa or the Higgs sectors. Within this class, we shall 
consider the following sub-cases: 

(a) The 2HDM-Va has Hermitian Yukawa matrices, but the Higgs sector is CP violating. To work out a concrete 
example we shall also consider a four-texture for the Yukawa matrices. 

(b) For the 2HDM-Vb we will not assume Hermiticity in the Yukawa matrices, while the Higgs sector is CP conserving. 
Again, to work out an specific example we shall consider the four-texture case for the Yukawa matrices. 



B. Solutions to the FCNC problem 

When both Higgs doublets couple to up- and down- type fermions, FCNSI are allowed An acceptable suppres- 
sion for FCNSI can be achieved with the following mechanisms: 
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Model type 


1 st family 


2 nd family 


3 rd family 


I 


Hi 


Hi 


Hi 


II 


H2 — > U quarks 
Hi —5- d quarks and leptons 


Hi — > u quarks 
Hi — > d quarks and leptons 


Hz — ¥ u quarks 
H2 —¥ u quarks and leptons 


III 


Hi and H2 


Hi and H2 


Hi and H2 


IV-*, b, r 


Hi 


Hi 


H 2 


V 


Hi and H2 


Hi and H2 


Hi and H2 



TABLE III: Higgs couplings to the fermion families. 



Universal Yukawa textures. Suppression for FCNC can be achieved when a certain form of the Yukawa matrices 
that reproduce the observed fermion masses and mixing angles is implemented in the model. This could be 
done either by implementing the Frogart- Nielsen mechanism to generate the fermion mass hierarchies [51( , or by 



studying a certain ansatz for the fermion mass matrices [37[. The first proposal for the Higgs boson couplings 
[3a ]. the so called Cheng-Sher ansazt, was based on the Fritzsch six-texture form of the mass matrices, namely: 

/ c q 

C* B q I ( 1 ) 

V B* A q 



Mi 



Then, by assuming that each Yukawa matrix Y± 2 has the same hierarchy, one finds: A q ~ m q3 , B q ~ ^/m g2 m q3 
and C q ~ ^/m qi m q2 . Then, the fermion-fermion'-Higgs boson (ff'°) couplings obey the following pattern: 
Hfifj ~ ^/mJ7nTJ~/mw, which is also known as the Cheng-Sher ansatz. This brings under control the FCNC 
problem, and it has been extensively studied in the literature to search for flavor-violating signals in the Higgs 
sector 

In our previous work we considered in detail the case of universal four-texture Yukawa matrices plj | , and derived 
the scalar- fermion interactions, showing that it was possible to satisfy current constraints from LFV and FCNC 
[52l [53j . Predictions for Higgs phenomenology at the LHC was also studied in ref . [54| HH ■ We can consider 
this a universal model, in the sense that it was assumed that each Yukawa matrix Y{ 2 has the same hierarchy. 

Radiative Suppression of FCNC. One could keep FCNC under control if there exists a hierarchy between Y"' 
and Y 2 ' d . Namely, a given set of Yukawa matrices is present at tree-level, but the other ones arise only as a 
radiative effect. This occurs for instance in the MSSM, where the type-II 2HDM structure is not protected by 
any symmetry, and is transformed into a type-Ill 2HDM, through the loop effects of sfermions and gauginos. 
Namely, the Yukawa couplings that are already present at tree-level in the MSSM (Y d ,Y 2 ) receive radiative 
corrections, while the terms (Y 2 d , Y™) are induced at one-loop level. 

• Alignment of the Yukawa matrices Another solution to the FCNC problem that have been discussed recently 
assumes that the Yukawa matrices could be aligned However, it seems that if such assumption holds 

at a high energy scale (much above the EW scale), it no longer holds at a low-energy scale [Hj|. 



III. THE LAGRANGIAN FOR THE 2HDM 



The most general structure of the Yukawa lagrangian for the quark fields, can be written as follows: 



quarks = fiyD^ + jfiyD ^jO + ^0 ytf^O + ? ytf^O + ^ (2) 

where Y^ 2 D are the 3x3 Yukawa matrices, denotes the left handed quarks doublets and ur, df> correspond to 

the right handed singlets. Here (£12 = io~2<t'i 2- The superscript zero means that the quarks are weak eigenstates. 
After getting a correct SSB [59l— f62j] . the Higgs doublets are decomposed as follows: 



"i+yi+»xi 

V2 



(3) 
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<^2 - y e^ 2 -H| 2 +i X2 J ■ ( 4 ) 

where the v.e.v.'s Vi and i>2 are real and positive, while the phase £ introduces spontaneous CP violation. Now, 
we transform the quarks to the mass eigenstate basis through the rotations: ul,r — Ul,ru^ r , <1l.r = Dl ro^ r , to 
obtain: 

quarks = u LUL Y»p+DyR + d L D L Yf?±±^D R dR 
+u L U L Y 2 D <p+D R d R + d L D L Y 2 D V2 ^1 X2 D R d R 
+u L U L Y^ j Lpl ^l Xl U R u R - d L D{Y^^U R UR 
+u L U L Y 2 u iP2 ^l X2 URUR - d L D[Y 2 u V2 U R u R 

+u L M u u R + d L M D d R + h.c, (5) 
Then, the (diagonal) mass matrices are given as follows: 

l^ = 2LY? + e-*^Y° (6) 

and 

M D = ^Y° + e*^% D , (7) 

where Y^ 2 = U L Y^ 2 U R and Y X D 2 = D L Y° 2 D R . Then, one can split the Yukawa couplings into the neutral and 
charged terms, both for the up and down sector. The neutral couplings for the up sector are given in terms of 
(four-components) Dirac spinors as follows: 

Cr™' = uY^^^P R u + uY 2 u ^^P R u 
v 2 v 2 

+uY^ fl±^P LU + uY 2 ^^±^P L u 
v 2 v 2 

+uM u u, (8) 

where Pl the chiral operators. In order to arrive to the final form of the Yukawa lagrangian we need to 

include the Higgs mass eigenstates. When one allows for the possibility of having CP violation in the Higgs potential, 
the CP even and CP-odd components get mixed [65|. This CPV Higgs mixing is included as follows, 







(*\ 




= R 




Xi 













(9) 



with H4 is a Goldstone boson. The matrix R can be obtain when by relating equations © and ^ with the physical 
Higgs mass eigenstate 

( G+v a +H+w a \ 



r 


qrl 




1 


COS #12 COS 6*13 


— sin 0i2 — % cos 012 sin #13 


2 


sin #i2 cos 6<i3 


cos #12 — i sin 012 sin #13 


3 


sin 6»i3 


i COS 013 


4 


i 






TABLE IV: Mixing angles for Higgs bosons which consider spontaneous and explicit CPV [63( ] . 
where a = 1,2, r=l,...,4 and z) Q , u> Q are the components of the orthogonal eigenvectors of unit norm 1 

v = ^ vi , V2 j = ( cos /3, sin /3 ) 



and 



fu>i, $2 J = ( — e ^sin/3, cos/3^ 



(11) 



(12) 



The values of q ra are written as combination of the Oij , which are the mixing angles appearing in the rotation matrix 
that diagonalize the mass matrix for neutral Higgs; table IIVI shows the different values for the q r 's. 



It is convenient to write the following relation, for a = 1, 2 

<pi + ixi = (** cos ? ~ ?r 2 e-^ 23+ « sin/?) ff r (13) 

r 

and 

^2 + iX2 = («ri e * sin ^ + 1r2e- lf>23 cos/3) i? r . (14) 

r 

Then, we arrive to the final general form of the neutral Higgs boson couplings for the up- type quarks: 



£uT ral = * ( S ijr + 7 5 ^",) UjHr + UiMV Uj , 



(15) 



with 



and 



q u 

ijr 



pu 

ijr 



^M^^-^-tan^^e^ 2 



1 



2\/2cos^ 

Similarly, for the down-type quarks we find: 



1k2 e Y 2ij 1k2e I 2i j 



)■ 



(16) 



(17) 



(18) 



1 Here we shall follow closely the notation of Haber and O'Neil [63j ] . 
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with 



ijr 



= -1m 



2v 



<3 



Qki + Qki ~ tan 



2\/2cos/3 



and 



1 



g fel - q* kl + tan/? (g^^ - fee^*^ 



9fc2e *2 9fc2 e X 2 



2y/2 COS /3 

On the other hand, the Yukawa couplings for charged states are given by 

„5 



-<Pi 



+h.c. 



+ I-7! ( ^ M U_ e i Han pyvA V _^ I _jl Y Uty 



(19) 



(20) 



(21) 



where V denotes the CKM matrix. The physical eigenstates for the charged Higgs boson (H + ) can be obtain 
through the following rotation: 



cos /3 -e^sin/3 \ ( G ± 



<P2 I \ e ±l ^ sin (i cos (i 
Therefore, the Yukawa couplings for charged Higgs are 



C 



H- 



cos/3 
-h.c. 



V2 

1 H+fvY? 1 -^ 



V2 



(22) 



(23) 



IV. SOME LIMITING CASES 

A. The THDM-V with explicit CP violation (2HDM-Va) 

In this case we assume the hermiticity condition for the Yukawa matrices, but the Higgs sector could be CP violating. 
For simplicity we shall consider that the Yukawa matrices obey a four-texture form, and CP is violated explicitly in 
the Higgs sector. 

As it is discussed in the appendix [K\ the assumption of universal 4-textures for the Yukawa matrices, allows to 
express one Yukawa matrix in terms of the quark masses, and parametrize the FCNSI in terms of the unknown 
coefficients Xij > namely K^j = Xij ^ mi J a3 , where the hermiticity condition reads Xij — x\j ■ These parameters can be 
constrained by considering all types of low energy FCNC transitions. Although these constraints are quite strong for 
transitions involving the first and second families, as well as for the b-quark, it turns out that they are rather mild 
for the top quark. 

Then, from (|16l) and (IT7|) . one obtains within 2HDM-Va, the following expressions for the couplings of the neutral 
Higgs bosons with up-type quarks, namely: 



1 jj jTYliVtlj 

Si ]r = 7T M *] Wri + Iri - tan^ (q r2 + q r2 )] + —7= -Xij (q*2 + Qri) 

2v 2V2wcos/3 



(24) 
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and 

P ior = ttMV [q* rl - q rl - tan (q* r2 - q r2 )\ + V - 3 X ij (£, - ft*) > (25) 

2v2ucosp 

similar expressions can be obtained for the down-type quarks and leptons, as well as for the charged Higgs couplings. 



B. The Yukawa Lagragian for the 2HDM-Vb 

In this case we shall consider that the Higgs sector is CP conserving, while the Yukawa matrices could be non- 
hermitian. Then, without loss of generality, we can assume that H3 is CP odd, while Hi and Hi are CP even. Then: 
cos 0i2 = sin (0 — a), sin 6*12 = cos (0 — a), sin #13 = 0, and e - ' 6 * 13 = 1. The mixing angles a and that appear in 
the neutral Higgs mixing, corresponds to the standard notation. The expressions (I10p for the neutral Higgs masses 
eigenstates can be written now in terms of the angles a and 0: 

$° = (v + h°sm(0-a) + H°cos{0~a) +iG°)d a 

+ -^= (h° cos (0-a)- H° sin (0 - a) + iA°) w a , (26) 

where a = 1,2, and for the CP-conserving limit v a and w a have a vanishing phase £ = 0. Additionally, when one 
assumes a 4-texture for the Yukawa matrices, the Higgs-fermion couplings further simplify as Y 2 ij = Xij"^j — • Then, 
the corresponding coefficient equation for up sector and h (r = 1) are 

= -Ml [sin(/3 - a) - tan/3cos(/3 - a)} + ^^Z^fc^l ^ + x f.) (27 ) 



ijl 

For H° (r = 2) one finds: 



p u _ V^hmJ cos{0 - a) ^ ^ 



2y/2v cos/3 



(Xij ~ Xlj) (28) 



S% 2 = [cos(/3 - a) - tan/3 sin(/3 - a)} + ^^ Sm ^ a) ( Xy + x j.) , (29) 



pu 

^ij2 - 



Vm^m- sin(/3 - a) , f x 
2^ cos/3 V X ^ ^ 



Finally, for A (r = 3) one obtains: 



*-««?(*-*)• (31) 



^"3 = tan/? - i-^ >— ( Xij + X U (32) 

2w J 2V2t;cos/3 V J ' 



C. The 2HDM of type I, II and III 

It is interesting, and illustrative, to consider the limit when the general model becomes the 2HDM-III, within 
2HDM-III, one has that the Yukawa matrices obey a 4-texture form, and also Yf — Yj , namely: 



Y 2ij - Xij - 
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Model type 


Hi 


H 2 


H 3 


2HDM-I 
2HDM-II 
2HDM-III 
2HDM-V 


s ^ 0, P = 
5 / 0, P = 
5 / 0, P = 
5, P^O 


5 ^ 0, P = 
S^O, P = 
5/0, P = 
5, P^O 


5 = 0, P/0 
5 = 0, P / 
5 = 0, P/0 
5, P^O 



TABLE V: Yukawa couplings for neutral Higgs. The 5 and P can be for up and down sector. 



r 


QU 


pu 

J -ijr 


O ^ j r 


pd 


1 

2 
3 


v sin p L J 

aina R M% 

v sin p L J 







i cot /3 7i^"£^ 


v sin p *J 

sina «M° 

v sin p *J 







£ cot p" n#-D 



TABLE VI: Explicit values of the Yukawa couplings for neutral Higgs in 2HDM-I. 



The condition of Hermiticity means then Xij = Xij ■ Within 2HDM III, we shall consider that the Higgs sector is 
CP conserving. Therefore the Yukawa couplings take the following form. For h° one gets, 

1 



5 1 " = -M%(sin(ft-a) + ta,nftcos(ft-a)) 

J v 

Xijy/rriimj cos(ft - a) 



(33) 



y/2v cos ft 

while 

P?ji = 0- (34) 

Then for H°, 



S"; 



ij2 



v 10 cos ft 



Xij^rriimj sm(ft - a) 



V2v 



cos ft 



(35) 



and 



P%2 = 0- (36) 

Then for A , 



QU _ r, 
°ij3 — u i 



(37) 



and 



r ij3 - 



I 

V2u cos (3 



(38) 



One can also reduce the general model to the 2HDM types I and II, by eliminating some of the Yukawa matrices 
Y^' D = and Y? = Y 2 D = 0, accordingly. The tables |V] and IVII summarize the corresponding results, they include 
the expressions for the neutral Higgs couplings with up and down type quarks, and similar results hold for the leptons. 
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r 


QU 

Oij r 


pu 

r ijr 


Qd 


pd 

r ijr 


1 

2 
3 


COB "MV 

v sin p "J 

sin a a mE 

v sin p "J 







1! L J 


sin ™ MR 

v cos p J 

cos a „ MR 

v cos p L J 







1) L J 



TABLE VII: Explicit values of the Yukawa couplings for neutral Higgs in 2HDM-II. 




h° 
P 




FIG. 1: Tree level Feynman diagrams for the decay. Right diagram is for h — > W be while left diagram is for h 



W~bc. 



V. PROBING THE CP VIOLATING HIGGS COUPLINGS THROUGH THE DECAY h 



cbW 



In this section we shall evaluate the asymmetry coefficient for the decay h — ¥ cbW in order to analyze presence of 
both FCNSI and CPV within the 2HDM-X. In the SM the FCNC are suppressed, but in the 2HDM extensions these 
processes are found even at tree level. We consider the neutral Higgs boson decay h — > Wbc at tree level. Two 
diagrams contribute to this decay, the first one is through the FCNC coupling h — > t c — > W~bc, its Feynman 
diagram is shown on left figure [T] The other one is through h — > W + *W~ — > W~bc, also shown in figure [1] 

The couplings of the neutral Higgs with the quarks and the charged boson W with the neutral Higgs are written 
as i (8231 +7 5 -P;>3i) an d igMy/qiig^" ', respectively. The other vertices are the usual SM contribution. The average 
amplitude for these diagrams is thus 



\M\ =\Mi 



■M\M 2 + M\Mi 



(39) 



We can obtain an approximation when the terms proportional to the charm and bottom masses are neglected. 
Then, the expressions for the squared amplitudes are 



\Mi 



4M^ 



\Pt(q)\ [4 IS&i " P 2 u 3 il Pi ' P2P1 ■ QP3 ■ q + 2 |5 2 " 31 - F 2 " 31 | 2 M^p 2 ■ qp 3 ■ q 



+ (V23I + P 23l\ 2 ™l ~ 1^231 ~ ^31 1' ? 2 ) (2P1 ' P2P1 • P3 + M^p 2 ■ p 3 ) 



(40) 



\M 2 \ = g 4 (qnY \V cb \ 2 \P w (k)\ 2 (M^ P2 ■ p 3 + 2pi • P2Pl ■ P3 ) , 



(41) 



M\M a = ^(Sr 



M w 



231 



PS 1 )q n V cb P*(q)P w (k) 



(M^p 2 



p 3 + 2pi ■ p 2 pi ■ p 3 ) 



(42) 



and 



M\Mi 



g i m t 



Mi 



w 



(Sj8l + ^23l) quV cb Pjr(k)P t (q) (M^p 2 ■ p 3 + 2 Pl ■ p 2Pl ■ p 3 ) 



(43) 



11 



where the W boson propagator is written in the Feynman-t'Hooft gauge Pw(k) — (k 2 — + iMw^w) 1 and 

P t (q) = [q 2 — m 2 + im t T t ) 1 . In order to find the asymmetry coefficient we also need to calculate the conjugate 
decay, that is, h — > W + bc. We denote the average amplitude as 



\M\ =\Mi\ + \M 2 \ +Mi M 2 + M 2 Ml (44) 
The square terms are the same as the above, 



Mi 



\Mi,2\ , while for the interference terms we have 



4 

M\M 2 = ^r 1 (5 2 " 31 + P 2 ll 31 ) q u V cb P t (q)P^(k) (M^p 2 ■ p 3 + 2p x ■ p 2P i • p 3 ) 



and 



„4, 



M\M X = ^ (S& + PftJ qil V cb P;(q)P w (k) (M^p 2 ■ p 3 + 2pi • p 2Pl ■ p 3 ) . (45) 
Then, the width for the decay is 

Th-^WSc = ^3 J J r (W* + + M\M 2 + MjM?) dxdy, (46) 

where the dimensionless variables are defined as x = and y = 2 -^ 2 -. All details about the decay kinematic were 

rn h y m h J 

included in the appendix |B] The definition for the asymmetry coefficient is 

A _ >W+bc ~ >W~bc jA7 , 

ACPV - r= — ■ (47) 

1 h — >W+bc "+" 1 h — >W-bc 

The final result, for the decay asymmetry is given by: 

2F c6 Re (q kl )Im + P? jk ) (J 10 + J12) 
A CPV (S% k ,P% k , Re (q kl ) , m h ) = ^ 3 -L , , (48) 



where 

/ (Sijki Pijk^ e (Qkl) 1 m h) = \Sfj k - Pi jk j ( Ji + J 2 - J4 - >h) + \Si jk + P"j k j ( J 3 + J 5 ) 

+<?Re {q kl f \V cb \ 2 {J 7 + J 8 ) + 2Re (q kl ) Re (Sf jk + P% k ) V cb (J 9 + J u ) . (49) 

The J's are integrals obtained from the decay kinematic, which are shown in appendix [B] as well as the others 
parameters defined in previous sections. 



A. Asymmetry in 2HDM-Va 

Let us discuss now the resulting expression for Acpv for two subcases within 2HDM-V. We fix i = 2 and j = 3 in 
equations (fM|) and (f25j) in order to obtain the appropriate parameters within 2HDM-Va, then we find: 



5 23i = 7T7k ^23 (9l2 + 912) (50) 

2v2ucos/i) 



and 



y/m c m t 
2\/2ucos/3' 



*&x = W V ' ° \ X2S (QI2 ~ Q12) ■ (51) 



Then, the asymmetry coefficient is 
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V / m c m t V c bX23 (Jw + J12) cos 2 6x2 cos 6x3 sin 13 

A 2 HDM-Va - r , (52) 

V2M w f {p, 612, 6»i3, X23, rn h ) cos p 



where 



2 

/ (/?, 012, 013, X23, m h ) = 32 ^r t C Q 2 s 2 ( Sin2 + C0S2 Sm 2 M (^1 + ^2 + J 3 ~ J 4 + J 5 - JJj) 

+ |V C& | 2 cos 2 12 cos 2 0x3 (J 7 + J 8 ) - V ^ T " tl/cbX2 , 3 cos0 12 cos0 13 sin0 12 ( J 9 + J u ) . (53) 

V2Afiv cosp 



B. Asymmetry in 2HDM-Vb 

The appendix lAl shows the four-texture structure for Yukawa matrices, nevertheless for practical evaluation of the 
asymmetry we write the texture parameter in Euler complex form as X23 — IX23 |e ,l/23 . Then, we fix the equations 
(|27p and (|28p for i — 2 and j — 3 in order to obtain the required element for 2HDM-Vb, 

cos(/3 - a)ym c m t / f \ 

*231 - o (X23 + X23 ) ( 54 ) 

2v 2v cos p v ' 

and 

_ cos(/3 - «)ym c m t / f \ 

^231 - 7T7k 2 l* 23 ~ ^23 J ■ ( 55 J 

2v2wcosp v ' 

Then, for this case the asymmetry coefficient is given by: 

gV cb ^/m c m t Ixas] sin ^23 cos(£ - a) sin (f) - a) 

A 2H DM-Vb 7=— — — (Jio + J12) , (5b) 

y/2M w cos Pf (a, P, X23,m h ) 



where 



tin \ gm c m t cos (p - a) 2 

/ (a,p,X23,w /l ) = 2 2^ — |X23 1 C-n + J 2- J4- J6 + J3 + J5J 

OiLilVlyy COS fj 

gy/m c m t Vcb , , cos ^23 sin (/} - a) cos(/3 - a) 

+ V2M ro 1X231 c^ (J9 + Jll) 

+ ff sin 2 (/3-a)|14 6 | 2 (J 7 + Ja). (57) 



C. Numerical results 



We shall discuss in detail the result for 2HDM-Vb. The asymmetry depends of the five free parameters. One of 
them is the Higgs boson mass which appears in the J integrals, the other ones are the mixing angles a, p and the 
complex parameter X23 • The mixing angle P is taken within the values 1 < tan /3 < 50 [64| . For the mixing angle a 
we study three scenarios, a < p, a ~ p and a > p. The phase V23 is fixed to the value 0.1 in order to analyzed a 
similar value to the phase from the CKM matrix. For each scenario we take two possible values for |x23|- Therefore, 
the scenarios studied here are: 

i) a < P for |X23 1 = 0.9 and |xas| =0.1, figured 

ii) a w P for |x2 3 | = 0.9 and |x2 3 | =0.1, figure [31 

iii) a > P for \x2s\ = 0.9 and |x2 3 | = 0.1, figure[U 

We use the reported values m t = 171.2 GeV, m b = 4.2 GeV, m c = 1.27 GeV, M w = 80.39 GeV and sxxx6 w = 0.231 

From figs [51 [3] and 2] we obtain asymmetry values of the order 10 3 to 10 2 (10 4 to 10 3 and 10 3 ) for scenario i) ( 
ii) and iii))within 2HDM-Vb. We have also analyzed the numerical results for the CP asummetry for case 2HDM-Va. 
We also find that the size of this asymmetry depens strongly on the phases. 
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FIG. 2: The asymmetry as function of tan/3 for scenario 1, on left for |x23 1 = 0.1 and on right for |x23 1 
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FIG. 3: The asymmetry as function of tan /3 for scenario 2, on left for |x23 1 = 0.1 and on right for |x23 1 = 
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FIG. 4: The asymmetry as function of tan/3 for scenario 3, on left for |x23 1 = 0.1 and on right for |x23 1 
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VI. CONCLUSIONS 



In this paper we have present a broad discussion of the most general formulation of the Two-Higgs doublet extension 
of the SM, which we name as 2HDM-X. Then, we have defined in a model named 2HDM-V, which has the possibility 
of including both FCNC and CPV, and have presented the corresponding Lagrangian for both the neutral and charged 
Higgs sectors. 

The limits when 2HDM-X reduces to one of the known versions (2HDM-I, II, III) has also been discussed; in these 
cases each pattern of Higgs- Yukawa couplings holds for all families. To identify the class of family non- universal 
models, we have used the label 2HDM-IV, where we include models where one Higgs doublet couples only to a certain 
type of fermion, for instance to the top quark or the third family, or to neutrinos only. 

Finally, we have also evaluated the CPV asymmetry for the decay h — > cbW, which allows to test the presence of 
both FCNC and CPV that associated with model V. We found that for certain optimal range of parameters the decay 
asymmetry could be of O(10~ 2 ) to O(10~ 4 . These asymmetry values for three scenarios were obtained in the case 
of the 2HDM-Vb. Similar results arise within 2HDM-Va. The asymmetry behavior has a dependency proportional 
to the mixing complex parameter X23- The mixing angles a and j3 control the shape of the graphs. The asymmetry 
keeps same shape for the Higgs boson mass range between 115 GeV and 160 GeV. 

In order to detect this asymmetry we could have to resort to a linear collider, since the final state seems difficult 
to reconstruct at a hadron collider. Although a final conclusion would require a detailed simulation study, which we 
plan to address in a future publication [65j. 
Acknowledgments . 
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Appendix A: 2HDM-III with four- Textures 

Here we shall summarize the result for 2HDM-III, namely we assume that both Yukawa matrices and Y% have 
the four-texture form and are Hermitic; following the conventions of [2o| . the quark mass matrix is then written as: 



M q = 



( ° 

V 



B„ 





, B q 
B* A n 



q —r one recovers the six-texture form. We also consider the hierarchy: 



when B 

1 4? I » I B i 

Because of the hermicity condition, both B q and A q are real parameters, while the phases of C q and B q , <&B q ,c q i 
can be removed from the mass matrix M q by defining: M q = P^M q P q , where P q = diag[l, e l ^ c i , e^* B i+' I ' c g )] 5 and 
the mass matrix M q includes only the real parts of M q . The diagonalization of M q is then obtained by an orthogonal 
matrix O q , such that the diagonal mass matrix is: M q = O^M q O q . 

The lagrangian (2) can be expanded in terms of the mass-eigenstates for the neutral (h°, H°, A ) and charged Higgs 
bosons (H^). The interactions of the neutral Higgs bosons with the d-type and u-type are given by (u, u' = u,c,t. 
and d, d' — d,s,b.), 



B q |,| C q |, which is supported by the observed fermion masses. 
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c q 



g I m d 
2 \ m w 



+ 



9 ( m d 



2 \ rriw 



cos a 

a °dd 

COS p 

sin a 



, V2 sin(a - p) ( m w , - , 
H ~ I 1 (*2 )dd' 



cos/3 



>dd' 



g cos p \ m d 

y/2 cos(a — ft) ( mw 
g cos P \ m d 



d'H° 



(Y 2 d ) dd , 



>u0 



d'ti 



2 \mw 

g ( m u 
+ | — I - 

+ § f-^.« 

2 \miy 



2 \mw 



.g cos P \m d J 



7 5 d'A° 



sin a 
sin/3 

cos a 
sin/3 

cot /3<5. 



\[2 sin(a — P) ( mw 



g sin/ 



/ rrO 



^2 cos(a-/3) / mvi/\ , 
««u' : — ~ I I (J- 2 )i 



+ 



g sin, 
5 sin/3 V m„ 



u'h° 



mw 



7 u 



(Al) 



The first term, proportional to S qq > corresponds to the modification of the 2HDM-II over the SM result, while the term 
proportional to Y 2 denotes the new contribution from 2HDM-III. Thus, the ff'<j>° couplings respect CP-invariance, 
despite the fact that the Yukawa matrices include complex phases; this follows because of the Hermiticity conditions 
imposed on both Y q and Y 2 . 

The corrections to the quark flavor conserving (FC) and flavor violating (FV) couplings, depend on the rotated 
matrix: Y 2 = O q P q Y 2 P^O q . We will evaluate Y 2 assuming that Y q has a four-texture form, namely: 



yq 



( Cf \ 

C? B q B\ 



SI 1,1 B q |,| C q 



\ Bf A\ J 

The matrix that diagonalizes the real matrix M q with the four-texture form, is given by: 



(A2) 



O a 



( 



A q (\l-\\){\%-\\) 



Vq 



-Vq 



\Vq 



'iq \ (a 
Ma 

V A q (X 



a?(a?-a,) 



(A q — X^)(A q — X^) 



(A 



EMM 



I \«\l{A q -\ — \ 
V A g (A«-A«)(Ag-A« ) 
' A«(A«-A q ) 



A^(A,-A;)(Ag-A g ) 
A q {\l-\\)(\l-\l) 



A,(A«-Ai)(A^-A^) / 



where to 9 =| Af |, mi, =| A^ |, ml =| A| |, and ry g = X\jm\ (q = u,d). With m u — m", 
= mf , m s = mf, and = mjj. 
Then the rotated form Y 2 9 has the general form, 



m c — m 2 , and m t = m 3 ; 



T q P q Y q P\O q 

( (Y 2 q )n (Y 2 q ) 12 fc'Jis 

(Y q ) 21 (Y q ) 22 (Y q ) 23 

V (Y 2 q hi (Y 2 q ) 3 2 (Y q ) 33 



(A3) 



However, the full expressions for the resulting elements have a complicated form, as it can be appreciated, for 
instance, by looking at the element (F 2 9 )22, which is displayed here: 



(Y q ) 22 = Vq {C q *e^ +Cfe-^] 



(A q -X q ) m\m\ 



\ q 

A 2 



A q m\ 



B qAq 



2 q 

mi 



A - \ q 
+A q q 2 



x 2 q 



[B q * e l * B « +B q 2 e 



^ l (A q -X q )(ml-A q ) 
1 ml -X q 



(A4) 
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where we have taken the limits: \A q \, m|, m\ 3> m\. The free-parameters are: B\, B q , A\, A q . 

To derive a better suited approximation, we will consider the elements of the Yukawa matrix Y\ as having the same 
hierarchy as the full mass matrix, namely: 



d = ^.W™ (A5) 




B\ = b^(A q -Xl)(ml-A g ) (A6) 

B\ = bl{ml-A q + \l) (A7) 
A\ = a\A q . (A8) 

Then, in order to keep the same hierarchy for the elements of the mass matrix, we find that A q must fall within 
the interval (m| — m\) < A q < m\. Thus, we propose the following relation for A q : 



A q = ml(l- f3 q z q ), (A9) 



where z q = m\lm\ <C 1 and < fi q < 1. 

Then, we introduce the matrix x 9 as follows: 



(«)„ - 



m q m q 

— x q 



v 



m q m q 



v 



-4'"^ ( A1 °) 



which differs from the usual Cheng-Sher ansatz not only because of the appearance of the complex phases, but also 
in the form of the real parts x% = \\ 



>~9 

ij i 



Expanding in powers of z q , one finds that the elements of the matrix \ q have the following general expressions: 
Xii - K - (4*^ + c\e-^ )]n q + [a\ + b\- {bfe^*« + b\e~^ )]/?, 



X \ 2 = (c q e-^ b\) Vq [a q + b\- (fere***. + 6^"**^)]^ 



X? 3 = (4-ble-^) Vq ^P q 
Xh = 4v q + [4 + bl- (bfe^ + b\e~^ )\p q 



tt 3 = (ble-^-al)^ 
xls = 4 (All) 

While the diagonal elements Xu are re& l: we notice (Eqs. 14) the appearance of the phases in the off-diagonal 
elements, which are essentially unconstrained by present low-energy phenomena. As we will see next, these phases 
modify the pattern of flavor violation in the Higgs sector. For instance, while the Cheng-Sher ansatz predicts that 
the FCNC couplings (^2)13 an d {Y 9 )23 vanish when a\ — b q , in our case this is no longer valid for cos$s q 7^ 1. 
Furthermore the FCNC couplings satisfy several relations, such as: |x| 3 | = |xf 3 |, which simplifies the parameter 
analysis. Similar expressions can be obtained for the lepton sector. 

Appendix B: Decay kinematics for h — s- cbW 

For sake of simplicity we introduce the dimensionless scaled variables 



and 
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[x, y, zj = ( 



_ ( 2Bi 2E2 2E^ 



With this notation we can write the energy conservation as 

x + y + z = 2. 



(B2) 



(B3) 



In the Higgs rest frame, we just consider the contribution of the /ii, because /ii >> ^2,^3, and find the momentum 
expressions 



Pi ■ P2 = 


^(x + v + m-i), 


(B4) 


Pi-Pz = 


2 

^(i-y-Mi), 


(B5) 


P2-P3 = 




(B6) 


pi-q = 




-^■(x + y + m-l), 


(B7) 


P2 ■ q = 


m h / IS 


(B8) 


P3-q = 


^(2-a-»), 


(B9) 


k 2 = 


™1 (! + Ml - 


(BIO) 


q 2 = 


m^(af + 2/-l). 


(BH) 



Now, the functions |P t (g)| 2 , |P[y(9)| 2 j P? (q)Pw(k) and Pt(q)P^(k) with the dimensionless variables can be written 



as 



\P{q)Y 



m\{ x + y-\- (if + ^r 2 ' 



and 



P t *(q)P W (k) 



Pt(q)Pw(q) = 



™h (1 - x) + Mi7 

(x + y - 1 - n) (1 - x) + y/JIJIi'-fT + i [ V //7T (1 - x) - ^/JT[j (x + y-1- /J,)) 



ml 



(x + y—l — fx) + fiT 2 (1 - x) + M17 



(x + y - 1 - n) (1 - x) + y/flfli'yT - i [^77T (1 - x) - ^/JT[j (x + y-1- //)] 



ml 



(x + y-1- /i) 2 + ^r 2 (1 - x) 2 + ^i 7 2 



(B12) 
(B13) 
(B14) 

(B15) 



here fi = T 2 = £fc and j 2 = with T t sa 1.28 GeV and T w w 2.14 GeV are the SM full width decay for top 
quark and W boson, respectively [64j . The three body decay rate is given by the formula 



dT 



h >Wbc 



2m/, 



d 3 pi 



(2ttY 2E X 



Using the delta function to perform the p~s integral and setting the polar axis along pi, we have 



d^2 



(2ttY 2E 2 



d'pt 



(2iry 2E 3 



(2tt) 4 (5 4 (p-pi -P2-P3) ■ 



r 



h >Wbc 



Til + r22 + Tis 



(B16) 



(B17) 
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where we define 



r " = 25^3 / / l^ 1 ' dxdy > ( B18 ) 



r22 =2^/X (B19) 



and 



r X2 = 



f f (7W t 1 7W 2 +M^i)dxd 2/ , (B20) 

J J R, x y 



256tt 3 

with the R xy region is defined by 

~ (2 - x - V^ 2 " 4/ii) <y<^(2-a; + V^ 2 1 4/xi) (B21) 

and 

2 V / Mi < x < 1 +fi x . (B22) 



We can obtain the next results whether we write equations (|40l) . (|4Tj) , (|42|) and (1431) with dimensionless parameters 
and substitute in equations (|B18|) . (|B19j) and (|B20|) . 



_ 2 5 2 m /l 

1 11 — 

(16tt) 3 



1^231 " ^23l! 2 (Jl + J2 ~ J 4 - J 6 ) + 1^31 + A" 3 l| 2 (-ft + J B ) , (B23) 



and 



_ ff^ \V cb \ 2 m h 

22 ~ 512^3 [ 7 sJ ' ( ' 

ru = l^23i+-P 2 "3ilg 3 giiK b m fe + + ^ ^ + 

25o7r d 



1^231 + ^231 1 3 3 9H ( Jg fl _ Ji() CQS e+Jnsin6 _ J 12 cos 0) . (B26) 



256tt 3 

The Ji integrals, for i — 1, 12, are given by 



j^^ff i*+v-K-m*+v+n-m2-*-v) (B27) 

Mii ii^ (x + y - 1 - /j) 2 + /^r 2 

J 2 = / I {X + y -^- 1){2 2 - X - y) d X dy (B28) 
J J_r xb (x + y-1- fi) + /^r 2 

Jz = ±t f iX + y - t * 1 - l)il - y - fXl) dxdy, (B29) 
Mii A XH (x + y - 1 - /j) 2 + /^r 2 

■h = M / / 7 + ^ dxd Vi ( B31 ) 

J J B.xy (x + y - 1 - fi) 
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FIG. 5: Graphics for the integrals of the fu. 
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FIG. 6: Graphics for the integrals of the T22 and Ti2. 



J 6 = 



(x + y-l)(l-x + ^) 
fi„ (x + y—1 — /i) 2 + ^r 2 

(1 - X + Ml) 



dxdy, 



*-//. 



J 7 = Mi / / -2 -dxdy, 

J Jr xv (1-x) + Mi7 

(i - x + mi) [(x + y-i- M ) (l - x) + VTwrr] 



(x + y-l- m) 2 + Mr 2 (l-z) 2 + Mi7 2 



dxdy, 



I Ir 



(l-x + mi) [y^T (1 - x) - ygl7 (x + ?y - 1 - m)] 



{x + y-l- /i) 2 + Mr 2 (1 - xf + M i7 2 



dxdy, 



Jn 



/H ./ ./;,■ 



(x + y + mi - 1) (1 - y ~ Mi) [(x + y-l- m) (1 - a;) + VMMi7 r ] 



(a: + y - 1 - m) 2 + M r2 (1 - x) 2 + Mi7 : 



dxdy, 



(B32) 
(B33) 
(B34) 
(B35) 

(B36) 

(B37) 
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JV2 



y + /ii - 1) (1 - y - m) [y/JIT (1 - x) - ^/JTxl {x + y-1- //)] 



(x + y - 1 - n) 2 + itf 2 (1 - xf 



dxdy. 



(B38) 



The graphics for the Ji, i = 1, 12, are shown in the figures [5] and [6] 
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